Doping- and size-dependent suppression of tunneling in carbon nanotubes 
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We study the effect of doping in the suppression of tunneling observed in muIti-walled nanotubes, 
incorporating as well the influence of the finite dimensions of the system. A scaling approach allows 
us to encompass the different values of the critical exponent a measured for the tunneling density 
of states in carbon nanotubes. We predict that further reduction of a should be observed in multi- 
walled nanotubes with a sizeable amount of doping. In the case of nanotubes with a very large 
radius, we find a pronounced crossover between a high-energy regime with persistent quasiparticles 
and a low-energy regime with the properties of a one-dimensional conductor. 
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During recent years there has been much interest in the 
investigation of the electronic properties of carbon nan- 
otubes (CN) [1]. The reduced dimensionality of these 
systems leads to the appearance of unconventional ef- 
fects, such as a suppression of the tunneling conduc- 
tance at low energy scales. This has been interpreted 
as a signature of the so-called Luttinger liquid behavior 
[2,3], characterized by the absence of electron quasipar- 
ticles in the spectrum. Evidence of a power-law behav- 
ior in the tunneling density of states at low energies has 
been obtained from measurements in ropes [4], individ- 
ual single- walled nanotubes (SWNT) [5] and multi-walled 
nanotubes (MWNT) [6]. 

In the Luttinger liquid picture, critical exponents of 
observables like the density of states are not universal and 
depend on the interaction strength. In the case of carbon 
nanotubes, this refers to the ubiquitous Coulomb inter- 
action. However, a precise determination of the interac- 
tion strength is precluded by the fact that the electron- 
electron interaction is actually long-ranged, and it can be 
assimilated to a coupling constant, only after introducing 
a suitable infrared cutoff in the singular expression of the 
one-dimensional (ID) Coulomb potential [7,8]. 

Furthermore, the shape of CN becomes relevant in set- 
ting the strength of the electron correlations. In the case 
of MWNT, measurements of the conductance refer usu- 
ally to the outer layer, whose electronic properties are in- 
fluenced by the interaction with inner metallic cylinders 
[9]. Also, MWNT use to be significantly doped, what 
leads to the presence of a large number of subbands at 
the Fermi level [10]. The contribution of a large number 
of modes at low energies has then an appreciable impact 
in the enhancement of observables like the tunneling den- 
sity of states. 

The purpose of the present paper is to study the com- 
bined effect of the finite length of the electron system, 
on the one hand, and of the number of subbands at low 
energies, on the other hand, in the Luttinger liquid de- 
scription of CN. With this aim, we adopt a renormaliza- 
tion group (RG) approach, which is well-suited to obtain 
the energy dependence of quantities like the quasipar- 



ticle weight. In order to cure the infrared singularities 
arising from the long-range interaction, we implement a 
dimensional regularization of the theory, writing formally 
all bare quantities slightly away from dimension D — 1 
[11,12]. This formal artifact has a real physical meaning, 
since the deviation D — 1 yields a measure of the finite 
size of the system, with the limit Z? — > 1 corresponding 
to the case of increasingly long nanotubes. 

The Coulomb potential l/|r| can be represented in 
three spatial dimensions as the Fourier transform of the 
propagator 1/k^. If the interaction is projected onto one 
spatial dimension, by integrating for instance the modes 
in the two transverse dimensions, then the Fourier trans- 
form has the usual logarithmic dependence on the mo- 
mentum, Vc{k) — (l/27r) ln(fcc/fc). In this approach, kc 
is the memory that the system keeps of the finite size 
of the transverse dimensions. We choose instead to inte- 
grate formally a number 3 — D oi dimensions, so that the 
long-range potential gets the representation 
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where c(D) =T{{D-l)/2)/{2y/^f-^. In the limit £> ^ 
1, the expression of the potential in momentum space 
becomes 
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This shows that 1/{D — 1) corresponds actually to the 
logarithm of the high-energy cutoff kc. In the RG frame- 
work, where all quantities are scaled down to low mo- 
menta, the dimensionless quantity 1/{D — 1) has to be 
traded for the logarithm of the length of the system, mea- 
sured in units of the nanotube finite radius. 

The great advantage of dealing with the above repre- 
sentation is that the Coulomb interaction gives rise clas- 
sically to a scale invariant theory at any dimension above 
D = 1, which allows to proceed with the RG program. It 
is only at D = 1 that the mentioned logarithmic depen- 
dence on the momentum leads to an imperfect scaling 
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behavior, making the use of the RG approach not quite 
appropriate. Then, we write the hamiltonian for the lin- 
ear branches of CN in the form 
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where the ai matrices arc defined formally by {ai,aj} = 
26ij and p(p) are density operators made of the electron 
fields ^'Qcr(p), with a labeling the Fermi point and a the 
spin projection. The sum in Eq. (3) runs then over the 
usual four modes present in a single- walled nanotube, but 
it may include also the contribution from a large number 
of subbands in the case of a doped multi- walled nanotube. 

We focus on the scaling properties of the model as the 
cutoff A is lowered, when a large number N of subbands 
contribute to the electronic properties down to the Fermi 
level. Each subband is labeled with a different quantum 
number, which corresponds to the momentum in the di- 
mension around the nanotube. For this reason, the dom- 
inant processes are those where each scattered electron 
remains in the same linear branch. The main effect of 
the interaction is to dress the bare electron propagator 
with the polarization of the N different subbands given 

by 
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where 6(D) = ^mE±^0r^^ 

ization (4) is the analytic continuation of the known re- 
sult for two linear branches with opposite chirality, which 
we take away from £> = 1 in order to carry out a con- 
sistent regularization of the Coulomb interaction. After 
dressing the interaction with the polarization (4), the 
electron self-energy is given by the expression 
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The low-energy properties of the theory are investi- 
gated by taking the limit A ^ 0, where the self-energy S 
turns out to have terms linear in uJk and k that depend 
logarithmically on the cutoff. This is the signal that the 
scale of the electron wavefunction Z^/'^ and the Fermi ve- 
locity Vf are renormalized at low energies. The divergent 
contributions to the electron propagator read 
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where g = 2Nb{D)c{D)^ /vf, K{D) = ^'^(^%%)tl'^\ 

and / (D) = 2B^(D+i)/2r(£)/2)6(D) • 

The usual RG argument is that the renormalized prop- 
agator G must be a finite quantity, so that the divergent 
dependences on the cutoff A have to be reabsorbed in 
the scale of the wavefunction Z^/'^ and the Fermi veloc- 
ity VF [14]. Under a differential variation of A, Z^/"^ is 
renormalized according to the equation 
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The renormalization of vf can be translated into that of 

the effective coupling g = 2Nb{D)c{D)e'^ /vf, since the 
electron charge e is not renormalized in our model. The 
RG equation for g becomes 
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We are now in a position to study the influence of the 
long-range Coulomb interaction in the limit D 1. For 
this purpose, we start by considering the RG equation 
(8). The series showing in the r.h.s. can be summed up 
at D = 1, the flow equation taking then the following 
form in the neighborhood of that point: 
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ft follows from Eq. (9) that, aX D = \, there is formally 
a line of fixed-points covering all values of the interaction 
strength. However, we still have to take into account that 
the effective coupling g{D) becomes singular in the limit 
D ^ 1. We can shufffe this divergence into the initial 
value of the coupling, go{D), which turns out to have 
near D = 1 the asymptotic behavior 
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Then, by matching the behavior of go{D) with that of 
Eq. (9), we observe that the D — 1 factor in the r.h.s. 
of the RG equation is not completely canceled out in the 
limit D ^ 1 . This shows in a rigorous way that the 
bare ID long-range interaction is at a RG fixed-point for 
arbitrary values of the interaction strength. 

We consider next the RG equation (7) for the electron 
wavefunction scale in the limit D ^ 1. The series in the 
r.h.s. can be also summed up at = 1, with the result 
that the scaling equation in that limit reads 
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The function in the r.h.s. of Eq. (11) coincides with the 
anomalous electron dimension found in the exact solu- 
tion of the Luttinger model [15,16]. This ensures that 
the RG approach is a sensible way to obtain the low- 
energy properties of the model. We have to bear in mind 
that our description introduces however two important 
differences with respect to the usual treatment of ID in- 
teracting electrons. On the one hand, the slight deviation 
of the model from D = 1 allows to control the approach 
to the bare long-range interaction as the length of the 
system is increased. On the other hand, we have also in- 
corporated the effect of the number N of subbands that 
contribute at low energies, in order to account for the 
influence of doping in MWNT. 

With the RG approach we can face two different exper- 
imental conditions, depending on the magnitude of the 
typical energy scale involved in the measurement pro- 
cess, as compared to the spacing between subbands in 
the nanotube sample. When the latter is larger than the 
temperature or bias voltage applied to the sample, we 
are in a situation where the number of subbands can 
be taken as constant along the RG flow. Otherwise, for 
large enough temperature or bias voltage, the number of 
subbands that contribute at the scale of the high-energy 
cutoff is a decreasing variable as A — > 0. We deal sepa- 
rately with the two instances in what follows. 

RG approach with a constant number N of subbands. — 
In transport experiments, the typical scale of tempera- 
ture or bias voltage lies usually below the scale of the 
spacing between subbands. This has been so, even in 
the measurements made in MWNT, where the spacing 
corresponding to a typical diameter d « 17 nm is as 
small as fivp/d « 29 meV. In these conditions, the only 
subbands that contribute to the properties measured ex- 
perimentally are those crossing the Fermi level. 

We apply then Eqs. (10) and (11) to confront the 
experimental results on the tunneling density of states 
gathered from different nanotube samples. Starting with 
the measurements made in SWNT, we take a number 
of subbands iV = 2 in the equations and adjust the de- 
viation from D = 1 in accordance to the length of the 
experimental sample. Following the argument below Eq. 
(2), we use the correspondence 1/{D — 1) « ln(L/d), L 
being the nanotube length and d the nanotube diame- 
ter. A suitable choice corresponding to the experiments 
reported in Ref. [5] is L/d ~ 10'', which gives in turn 
D w 1.14. 

The measurements of the tunneling density of states 
in SWNT have shown a power-law dependence on en- 
ergy, with values of the critical exponent a accumulating 
around « 0.35 [4,5]. We have represented in Fig. 1 the 
estimates obtained from the r.h.s. of Eq. (11), for small 
values of _D — 1. The values of a have a smooth de- 
pendence on the length L of the system and fall around 
a w 0.35 iov N = 2 and D « 1.14, with a reasonable 
choice of the couphng c^/tt^wf ~ 1-5. 

Our results show also an overall agreement with the 
exponents measured in MWNT. It has been noticed that 



such systems arc significantly doped, so that a large num- 
ber of subbands are found at the Fermi level. The exper- 
imental conditions refer to a situation where « 5 — 10 
(without taking into account the spin degeneracy). It has 
been reported that the values of the critical exponent a 
measured in 11 different samples range from 0.24 to 0.37 
[6]. This variation can be accounted for within our RG 
approach by assuming that the number of subbands used 
in the renormalization may shift from A'^ = 2toA'^=10. 
Part of the drift observed in the critical exponent may be 
also due to the smaller aspect ratio of MWNT, although 
this fact is difficult to assess given the lack of information 
about the total length of the experimental samples. 
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FIG. 1. Estimates of the exponent a from the r.h.s. of Eq. 
(11), for eVw^VF « 1.5 . The different curves correspond, 
from top to bottom, to D = 1.14, 1.16 and 1.18 . 

We show then that the suppression of tunneling in 
MWNT can be softened by increasing the doping level. 
Our results may be also relevant in general for nanotubcs 
of large radius where there are a large number of sub- 
bands crossing the Fermi level. These instances can be 
considered as midway in the process of making contact 
with the physical properties of a graphene sheet. This 
requires taking systems with larger transverse size, what 
in turn may lead to a situation where the spacing be- 
tween subbands is smaller than the typical energy scale in 
the experimental measurements. One has then to change 
slightly the computational scheme, as discussed below. 

RG approach with a cutoff- dependent number of sub- 
bands N{A). — In samples of a very large diameter, one 
may envisage conditions where the temperature or the 
bias voltage are much larger than the subband spacing 
in the nanotube. The RG approach can be still imple- 
mented, but taking into account that the number of sub- 
bands contributing in the partial integration of modes at 
energy A depends on the value of the high-energy cut- 
off. Let us suppose for simplicity that the system is at 
half-filling, with the typical structure of two subbands 
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crossing at the Fermi level. For not too large energies, 
the number of subbands crossing the energy level e has 
then a linear dependence on e, -/V(e) = A'q + UpS, where 
A^o = 2 is the number of subbands at the Fermi level. 

We obtain the scaling of the quasiparticle weight Z 
and the different observables by imposing a dependence 
of N on the cutoff A according to the above formula. 
This description yields a sensible prediction for experi- 
ments where the average effect of a large number of sub- 
bands is measured. The dependence of the quasiparti- 
cle weight Z on energy as the cutoff is sent towards the 
Fermi level is represented in Fig. 2, for different values of 
Up. In the present instance, the dominant contribution to 
the power-law behavior of the tunneling density of states 
comes from the number of subbands varying with the en- 
ergy scale. This sets the value of a close to 1 from the 
start, what is the natural way of recovering the charac- 
teristic linear density of states of a graphite layer from 
our ID point of view. 
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FIG. 2. Energy dependence of the quasiparticle weight Z 
at dimensions 6/5 and 8/7, for different values of Up. 

We observe from the results in Fig. 2 that the quasi- 
particle weight Z tends to have a flat behavior at high 
energies for large values of the density of subbands Up. 
This is in contrast to the rapid decrease signaling the 
typical power-law behavior for small values of Hp . In the 
curves for Up = 100 and Up = 1000, we see the existence 
of a crossover between a regime with persistent quasipar- 
ticles and another characteristic of the Luttingcr liquid 
behavior. The physical interpretation is that, for high 
energies above the crossover scale, the system has simi- 
lar properties to the 2D graphene, while one has to look 
at sufficiently small energy scales (or large length scales) 
to measure the properties of the ID wire. 

We have obtained several results that may be checked 
against future measurements carried out in MWNT and 
nanotubes of very large radius. We predict that the ex- 



ponent giving the power-law behavior of the tunneling 
density of states may suffer a significant reduction upon 
doping those systems, with the possibility of reaching val- 
ues as small as 0.1 for « 40. We have also shown that, 
when dealing with nanotubes of very large radius, there 
is a high-energy regime with persistent electron quasipar- 
ticles which has properties closer to 2D graphene than to 
the Luttinger liquid. We believe such features may be of 
interest when developing carbon-based devices made of 
graphene and nanotube structures with different shapes. 

This work is partly supported by the Italian Research 
Ministry, National Interest Program. 
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